TPBU-11-95 



:ond-mat/9601083 
December 1995 

Upper and lower bounds on the partition 
function of the hofstadter model 

Alexander Moroz^ 

School of Physics and Space Research, University of Birmingham, Edgbaston, 

Birmingham B15 2TT, U. K. 
and 

The Erwin Schrddinger International Institute for Mathematical Physics, 
Pasteurgasse 6, A-1090 Wien, Austria 



ABSTRACT 

Using unitary equivalence of magnetic translation operators, explicit upper and lower 
convex bounds on the partition function of the Hofstadter model are given for any rational 
"flux" and any value of Bloch momenta. These bounds (i) generalize straightforwardly 
to the case of a general asymmetric hopping and to the case of hopping of the form 
tjniSj" + Sj~ n ) with n arbitrary integer larger than or equal 2, and (ii) allow to derive 
bounds on the derivatives of the partition function. 
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1 Introduction 



In what follows, we shall consider the Hofstadter model |1[ described by the Hamiltonian 

H = t 1 (S 1 + S* l )+t 2 (S 2 + S*), (1) 

where unitary operators Si and 5*2, sometimes called magnetic translation operators, 
satisfy the commutation relation 

S 1 S 2 = ooS 2 S 1} (2) 

with 

us = exp(— 2iiia), (3) 

a being a real number. In general we shall allow for additional hopping terms of the 
form tj n (S™ + Sj~ n ) with n arbitrary integer larger than or equal 2. However, as will 
become clear later, our results generalize straightforwardly the symmetric case t\ = t 2 
with nearest-neighbour hopping. Therefore, for the sake of clarity, we shall firstly derive 
our bounds for t\ = t 2 and then the relevant generalization for the case t\ ^ t 2 and 
additional hopping terms is discussed afterwords. 

The model ([!]) appears in many physical problems. Originally, it appeared in the 
study of the behaviour of electrons moving simultaneously in a periodic potential and an 
uniform magnetic field B M. Subsequently, it has been also used in early discussions of 
the integer quantum Hall effect 0, to describe superlattices and quasicrystals ||, and 
flux phases in the high-T c superconductors || . In the case of lattice electrons moving in 
an uniform magnetic field B, the parameter a(= &/&o) is the flux per unit cell measured 
in units of flux quantum $o(= hc/\e\), and, in the Landau gauge A = (—By, 0, 0), 

S t = e^ pl+ ~- By ) = e dl+2lviay , (4) 

and 

S 2 = e ^ 2 = e d2 . (5) 

In the case where a = p/q is rational with p and q relative prime integers, the model is 
known to enjoy SL(2, Z) symmetry ||, and, because of the periodicity over an enlarged 
unit cell f7|, one can define Bloch momenta, k = (k\, k 2 ), and a corresponding magnetic 



Brillouin zone ffl. In momentum space, the Hamiltonian H and operators Si and S2 can 
be written as q x q traceless matrices, 
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Despite its simple form and the extensive numerical knowledge available, the Hofs- 
tadter model defies exact solution (see Ref. || for a recent progress). Here, we shall 
concentrate on the partition function, 



Z(/3, ki,k 2 ) = Tre- mkl ' k2 \ 



(9) 



where (3 = 1/T and T is temperature. The total partition function is then given by the 
relation 

d 2 k 



(10) 



where the integral is over the Brillouin zone. The calculation of the partition function 
Z(/3, ki, k 2 ) can be viewed as a technical tool to calculate the density of states p(E, ki, k 2 ) 
H, because Z(f3, hi, k 2 ) is nothing but the Laplace transform of p(E, ki, k 2 ), 



POO 

Z(0, h, k 2 ) = / e-? E p(E, k^ k 2 ) dE. 
Jo 
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2 Results 



We are not yet able to calculate the partition function Z(/3, ki, k 2 ) exactly, however, 
we shall formulate upper and lower bounds on Z(f3,ki,k 2 ), which hold for any rational 
a — p/q. To derive the lower bound on Z(f3, ki, k 2 ), we shall use the Peierls variational 
principle JlO| which states that 

Z(J3) = Tr exp(-(3H) > £ e -^(*».»*n) j ( 12 ) 

n 

where {$ n } is an arbitrary set of linearly independent normalized vectors in a Hilbert 
space. Eventually, to obtain the upper bound on Z(/3, ki, k 2 ), we shall use the Golden- 
Thompson inequality ||11|| , 

Tre A+B <Tr(e A e B ), (13) 



valid for self-adjoint operators A and B. 

Abstract inequalities (|12|) and ([OJ) are well known, however, it is difficult to calculate 
the traces involved and to obtain explicit analytic bounds in general case. In the present 
case, by taking A = Si + S* and B = S 2 + S 2 , we shall calculate e A , e B , and the traces 
in Eqs. ( ]Ti2| ) and (|T^) explicitly and we shall show that for all /3, k±, and k 2 , 

q < max (Tr e~ pA , Tr e~ pB ) < Z([3, k u k 2 ) < - (Tr e~ pA ) (Tr e~ pB ) , (14) 



where 



and 



Tr e~ 0A = £ Q(h, I), Tr e^ B = £ Q(A&, I), (15) 



1=1 1=1 



Q(k,l) = exp{-2p cos[k + 2na(l - 1)]}. (16) 

Bounds (0), which constitute the main result of our paper, are valid for (3 of both sign, 
positive and negative, and they are saturated for (3 = 0, i.e., at infinite temperature. 
The upper bound on Z in ([14]) follows from (|13"D using the clustering property of Tr 
e A e B _ (i/g^Xre^Tr e B to be derived below. This clustering property is particular to the 
Hofstadter model and does not depend on f3 and the hopping parameters. 



Both upper and lower bounds on the partition function in (I4"D are convex functions in 
variable f3 (and hopping parameters t\ and t 2 ), as the partition function itself. Given con- 
vexity of e~ x , the first inequality (from the left) in (|T3|) follows from the Jensen inequality 



12] and the fact that J2i=i cos [k + 27ra(Z — 1)] = for all k. Using simple arguments 
presented in ]13[, bounds (|b|) are used in Sec. 4 to derive bounds on the derivatives of 
the partition function with respect to /3, ti, and t 2 . 



3 Calculation of traces and derivation of bounds 

The main step in calculating traces in Eqs. ( |i~2|) and (p~3|) and deriving bounds ([14]) is 
to use the unitary equivalence of operators Si and S 2 defined by relations (0) and (H). 
Indeed, one has 

S 2 = US X U~\ (17) 
where U is the unitary matrix with elements given by 

U lm = — a; . . . = J_ w -0-l)(m-l)+m(m-l)/2_ /^x 

By direct calculation, 

[t^tr 1 ]* = j:u u u* kl u-^ = lj; w H)M w -(i-i) = ^ (19) 

Z=l ^ i=l 

One can verify that Ui +qt7n = U^ m and 

£W 3 = I godd (20) 

' q \ ~Ui, m , q even. v ; 

For all j and / (summation over repeated indices suppressed), 

UflUji = ~ q ■ (21) 



In general, 



JJ*JJ. k _ l w -(3'-l)(*-() ^ ^l*-'! — V = - w _ U _l }(fe _ <) + S g nCfe _i )x | fe _ <K | fc _i| + l )/2 ^ 

The action of U on the Hamiltonian iJ, UHU^ 1 , interchanges momenta k\ and k 2 and 
implies that the spectrum of H is symmetric under this interchange. It also allows to 
show that the spectrum of H is equivalent to the spectrum of the Hamiltonian H of the 
form 

H = U-^DxU 1 ' 2 + U 1 / 2 D 2 U~ 1 / 2 , (23) 



lxsgn(fc-2) 1 



where D\ and D 2 are real diagonal matrices, (D)i m = 2cos(fc + 2iral)6i m , with subscript 
labelling the dependence on either k\ or k 2 . 

Due to (0), both Si and S 2 have the same spectrum £ [JTJ] which can be read off from 
(g), i.e., 

£ = {c^|j = 0,l,...,g-l}. (24) 

Normalized eigenvectors of Si are (<&j)i = Sji, and corresponding normalized eigenvectors 
of 62 are 

= — (l, uj- j , LU~ 2j , . . . , u-( q - 1)j ) . (25) 



y/Q 

One can verify that, in the standard scalar product in C q , ($j, $/) = \l/;) = (5^. The 
lower bound on in (|14]) then follows from the Peierls variational principle (12) by 



taking the complete set of vectors to be the set of eigenvectors of either 5*1 or S 2 . Indeed, 
one has 

H$j) = 2 cos(/ci + 2naj), (%, HVj) = 2 cos(£;2 + 2naj). (26) 
The same lower bound can be also obtained from the Bogoliubov variational principle 



- pF > In Tr exp(-/3H) - (3 {Tr [exp(-(3H)(H - H)] /Tr exp(-/3#)} , (27) 

where if is a trial Hamiltonian which is not necessarily a dynamical operator. One uses 
the Bogoliubov variational principle ( P7| ) with if to be chosen either of the operators A 
and B. In the case of general Bloch momenta k\ and k 2 , one has 

e A _ e Si+5* _ ^ ia g | e 2cos[fci+2 7 r(/-l)Q]| ^ g B _ gS 2 +S* _ f/^iag | g 2 cos[fc 2 +27r(«-l)a] | _ 

(28) 

The second term on the right side of Eq. fl2T|) then vanishes [jTBI and the remaining term 
gives the lower bound in (|T4| ) on the partition function Z(/3, ki, k 2 ), valid for any rational 
a and for all (3, k\, and k 2 . 

As for the upper bound flH) on Z(f3,kx,k 2 ), using Eqs. (ES) and (EI]) one has 



Tr(e-^e-^) = l)Q(k 2 , j)U*U tj = - £ 5(^1, QQ(k,j) 

= -(Tre^) (Tre-^), (29) 



which is the clustering property mentioned above. Combining Eq. ( p9|) with the Golden- 
Thompson inequality (|13|) then gives the upper bound ( |T4] ) on Z(f3, ki, k 2 ), 



Z(/3,h,k 2 ) < - (Tre-W) (Tre"^) , (30) 



q 

valid for any rational a = p/q and for all f3, k\, and k 2 . 

4 Further generalizations 

For a given (3, one can check numerically that [after dividing (|14]) by q] the upper and 
lower bounds on Z(f3, k\, k 2 ) depend only very weakly upon q and Bloch momenta. The 
dependence of the bounds on (3 is much stronger. As (3 increases, 

max (Tr e~ pA , Tr e' 1313 ) /Z(J3, h, k 2 ) « e~ 2/3 -> 0. (31) 

Therefore, the lower bound is far from being optimal and, in principle, it can be further 
improved using the Peierls variational principle ([T^) with a better choice of the complete 
set of vectors than used here. With respect to large (3 behaviour, the upper bound ( |I~4"D 
on Z(f3, ki, k 2 ) is much better because it mimicks the large (3 behaviour of Z((3, hi, k 2 ). 
The upper bound can also be further improved: by repeatedly using inequality 



Tr (XYf m+1 < Tr (x 2 Y 2 ) (32) 



for nonnegative matrices X and Y, with m integral and > 0, together with the Trotter 
formula, 

e A+B = lim (e A ' n e B l n ) n . (33) 

One finds for 2 < I < m that 

Tr (e-^e"^) > Tr (e^ A ' 2 e^ B ' 2 ) 2 > Tr ( e -/W* e -/^y 

> ( e -^/- e - /3B / m ) m >Tre- /3H = Z(/5). (34) 



Note that relation ( P2"D enables us to calculate explicitly 

uS gfl( m -{)x\m-l\(\m-l\+l)/2Q^ ki) l)Q(k 2i m -l), (35) 



-^(Si+S" 1 ) -P{S2+S^) 



I in 



where Q(k 2 ,m — I) is a discrete Fourier transform, 

Q(h, m-l) = -J2 ^- (J - I){m - Z) Q{h,j). (36) 

Obviously, in the asymmetric case, ti 7^ t 2? and even in the case where hopping has the 
form tjn^S™ + 5j~ n ) with n > 2 arbitrary integer, one can still use the unitary equivalence 
(|17|) with the matrix {7 given by Eq. (|T8|). One writes H = A + B, where matrix A (B) 
now includes all hopping terms in the x-direction (^/-direction) and proceeds with matrices 
A and B as before, i.e., uses them as the trial Hamiltonian H in the Bogoliubov variational 
principle ( $I7\) and in the Golden-Thompson inequality (|T3"D. Inclusion of hopping terms 
tj n (S™ + S~ n ) in the Hamiltonian H does not change the form of bounds ( |T4| ) and involves 
only modification of Q(k,l). For example, in the asymmetric case, t\ 7^ £2, with the 
nearest-neighbour hopping, the only change in bounds (14) consists in replacing A by t\A 
and B by t 2 B. 

5 Bounds on derivatives of Z(f3, ^1,^2) 

The partition function Z(/3, fci, /C2) as a function of j3 and the hopping parameters ti and £2 
is a convex function. Adapting simple arguments presented in [|13[ for a concave function, 
one can use bounds ( |T4"D to derive bounds on the derivatives of the partition function with 
respect to /3, ti, and £2- Indeed, let f{v) be a convex function, i.e., f"(v) > 0. Let us 
denote by / ± (f ) the upper/lower bound on / at the point v. Then 

^<»°> - LM < />o) < z+W-rw (37) 

v -vi v - v 2 

where Vi < v < v 2 . Inequality (^) then implies, for example, 



A^{!^ e ^ (Tre-^) -max(Tr e -^,Tre-^)| < ^(Po,h,k 2 ) 



< j^j j ~ (Tr e- poA ) (Tr e"* 5 ) - max (Tr e-* A , Tr e~^ B ) j , (38) 



where f3\ < (3q < (3 2 . 
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6 Conclusions 



To conclude, we have calculated traces involved in general bounds on the partition function 
[Eqs. ( |r2|) and (0)] and established explicit upper and lower bounds on the partition 
function of the Hofstadter model [Eq. flT4l)1 and its derivatives [Eq. Q59)], valid for any 
rational "flux" a, inverse temperature (3, and Bloch momenta k\ and k 2 . Bounds (0) are 
saturated for (3 = and, in the asymmetric case, if either t\ or t 2 approaches zero. 

Bounds ( |i4|) also imply constraints on the density of states of the Hofstadter model and 
supply constraints imposed by the coefficients of the secular equation for the Hofstadter 
model, 

det(AI - H) = X q + a q . 2 X q - 2 + a 9 _ 4 A 9 " 4 + . . . + det H = 0. (39) 

Here, depending on whether q is odd or even, either det H = 2[cos(g&i) + cos(g&;2)] or 
det H = 4 — 2[cos(qki) + cos(qk 2 )}. One can show that a q -(2i+i) = 0, a q ^ 2 = — 2 9, 

a q -4 = q{2q — 7) — 2q cos(27ra), 
2 

a q _ 6 = — q{2q 2 - 21q + 58) + 4q(q - 6) cos(2vra) - 4q cos(4vra), (40) 
3 



17| and that except for the constant term, det H, none of the coefficients aj depends on 
the Bloch momenta. 
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hospitality. This work was supported by the UK EPSRC Grant GR/J35214 and by the 
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